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>■ ' Abstract 

lO ' Lattice theories that contain chiral multiplets of fermions can have complex 

fermion determinants. This is for example the case for the U(1)l®U(1)r symmetric 
Yukawa model with mirror fermions, if the number of generations of fermions and 

(N ■ 

Q\ • mirror fermions is odd. Whether a numerical simulation of such a model is possible 



depends on the magnitude of fluctuations of the complex phase factor of the fermion 
determinant. We investigate the fermion determinant of the U(l) Yukawa model 
with mirror fermions for a physically relevant choice of parameters. The argument 
' of the complex phase turns out to fluctuate only very little and is at most of the 

■ order of 2 • 10~^. 

X 

^ 1 Introduction 

The problem of formulating theories with a chiral fermion content on a lattice has not 
yet been solved in a satisfactory way. As is well known lattice fermions are accompanied 
by unwanted doublers which one likes to remove at least in the continuum limit. Their 
appearance is related to the Nielsen-Ninomiya theorem which requires that to every 
lattice fermion with a given set of quantum numbers a mirror fermion with the same 
quantum numbers but opposite parity exists. For a naive discretisation of the fermion 
action half of the doublers play the role of mirror fermions 0. 

The Wilson term P], which achieves the removal of doublers in vectorlike theories, 
spoils chiral invariance and cannot be used directly in chiral theories. For fermions inter- 
acting with scalar fields via Yukawa interactions a chiral invariant Wilson- Yukawa term 
has been proposed [Q, but there is evidence that it does not fulfill its goal p. On the 
other hand it is possible to remove the doublers in chiral Yukawa models with the help of 
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a Wilson-like term but at the cost of introducing extra mirror fermions . This yields the 
minimal possible proliferation of extra fermions. If one does not like to have the mirror 
fermions in the physical spectrum, even at a higher mass scale, the task is to find out 
whether the mirror fermions can be removed or decoupled in the continuum limit. 

If gauge fields are neglected, chiral Yukawa models with mirror fermions are for a 
particular choice of parameters invariant under a Golterman-Petcher shift symmetry ^ 
with respect to the mirror fermion field p| . As a result the mirror fermions decouple from 
the fermions and from the scalar field. It remains a truly chiral set of fermions interacting 
with the scalar field. 

For a numerical simulation of such theories with the Hybrid Monte Carlo algorithm, 
however, the fermion fields have to be duplicated again. This is due to the following fact. 
The part of the action bilinear in the fermion fields can be written with the help of the 
"fermion matrix" Q{(f))yx as 



contains fermion fields ijj and mirror fermion fields Xi ^ind is the complex scalar field. 
The matrix Q is not necessary positive definite, while the hybrid Monte Carlo algorithm 
requires a positive definite fermion matrix. A duplication of fermion fields by introducing 
two generations '^^^^ and \E'^^-' with opposite chirality amounts to replacing Q by Q^Q 
and ensures the positive definiteness. This duplication is of a purely algorithmic origin 
and has nothing to do with the fermion doubling problem. Therefore it is desirable to 
find ways to simulate these models without additional duplication. 

In particular, for a U(1)l(S>U(1)r symmetric model the determinant of Q is in general 
complex, preventing a simulation of this model with standard Monte Carlo algorithms. 
Nevertheless a simulation without duplication would be possible depending on the fluc- 
tuations of the phase angle a of the fermion determinant. Two possible situations can be 
imagined: 

1. the angle a fluctuates only very little about zero. Then a simulation with | detQ| 
appears feasible and the phase factor expia can be put into the observables or 
neglected completely. 

2. a fluctuates strongly and a simulation without taking the phase of det Q into account 
is not possible. 

In order to find out which of these two possibilities holds we have investigated the fermion 
determinant in the U(1)l ® U(1)r symmetric Yukawa model with mirror fermions. 
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2 Calculation of the fermion determinant 



The U(1)l ® U(1)r symmetric Yukawa model with mirror fermions has been considered 
in ref. [^j in detail and we refer to this article for the definition of the action and the 
parameters. The 8 ® 8 matrix Q is given in a 2 2 block notation by 
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Here x and y are lattice points, the sum runs over eight directions of the neighbours, 
and fi is the unit vector in the direction of fi. The Euclidean 7-matrices are expressed in 
a chiral basis as 
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(4) 



and S4 = S4 = 1, where cri,2,3 denote the 



For yU = 1,2,3 we have = 
Pauli-matrices. For negative indices the definition is given by = — The Yukawa 
couplings of the fermions and mirror fermions are and G^, respectively. The term 
proportional to Kr is a chiral invariant Wilson-term which serves to give the fermion 
doublers masses of the order of the cutoff. The Wilson parameter r is set equal to 1 
usually. Finally K is the fermionic hopping parameter, whose critical value is 1/8 for 
vanishing Yukawa couplings. 

The hybrid Monte Carlo algorithm requires the inversion of Q, which is done by 
conjugate gradient or minimal residue algorithm. A calculation of the determinant of Q 
is more demanding. Straightforward application of standard determinant algorithms is 
not possible owing to the need of large storage. We proceeded in the following way. Let 
Li, z = 1, 2, 3, 4 be the size of the lattice in the direction i. We consider even Lj and set 
L4 = 2n. Q can be considered as a 2ri x 2n matrix consisting of block matrices belonging 
to time slices. Each block itself is a M x M matrix with M = 8L1L2L3. In a first step 
the matrix is reordered by an even number of exchanges of rows and columns such that a 
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band structure is obtained: 



mr(i) b 
a mT{2) 
b 
a 
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mT{3) 

b 



b 






b 
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mT(2n-i) a 
±a b mT{2n) j 



(5) 



The blocks a and b are constant, whereas the vrtj depend on the scalar field 0. The 
indicated signs near the lower right corner depend on the temporal boundary conditions. 
T(j) is a permutation of the time slices defined by 



T(l) = n, r(2n) = 2n, 
j, 7'(2j + l)=n + j forj 



T{2j)=n-j, T{2j + l)=n + j for j = 1, . . . , n - 1 . (6) 

The matrix above is now dealt with by LU-block factorisation |]T^. The resulting upper 
triangular matrix has blocks on its diagonal, which can be determined recursively from 
the blocks rrij, a and b. The complete determinant is given by 



2n 



det Q = Y[ det . 



(7) 



k=l 



The "little" matrices are not sparse. But for a ■ lattice they are small enough 
such that their determinants can be calculated by standard numerical routines. Similar 



techniques have been used by Toussaint as indicated in [[11 



The storage needed does not depend any longer on L4 and amounts to 6144(LiL2-L3)^ 
bytes using complex extended precision. In our case this is slightly more than 3 Mega- 
words, which means a reduction by a factor of about 1/10 compared to a standard treat- 
ment. The CPU time is proportional to L4 instead of which standard routines require. 
For L4 = 8 it takes 3 minutes to calculate one determinant on the Cray Y-MP with a 
high degree of vectorisation (300 MFlops). 

Because the correctness and precision of the algorithm is crucial we have tested it in 
different ways: 

1. For lattices smaller than 4*^ a comparison with the results obtained with the help of 
IMSL-library routines has been made. 

2. For constant scalar fields the determinant can be calculated exactly by Fourier 
transformation, and we compared with the resulting values. 
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The relation 



A. 
dX 



[lndet(g + A7„„)]A=o = <5; 
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with {Iuv)xy = Kx^vy ^nd variable scalar field was verified numerically, where the 
left hand side was evaluated by our determinant algorithm and the right hand side 
with the help of the conjugate gradient algorithm. 

4. For particular choices of the Yukawa couplings the determinant has to be real. We 
checked this property to a high precision. 

In all cases the deviations were of the size of the numerical precision of the computer. 
For example in the last mentioned item the ratio of the imaginary to the real part of the 
determinant was of the order of 10~^^. 

3 Results on the fermion determinant 

Some general results on the fermion determinant in the U(1)l<8)U(1)r symmetric Yukawa 
model are available. 

1. The expectation value detQ is real. This is due to dct Q{(f)'^) = (det(5(0))*, and 
the weight of a scalar field (j) in the functional integral being the same as that of its 
complex conjugate 0"*". 

2. For G-^ — the determinant is real for each single scalar field configuration, as 
can be shown using the symmetries of the action. 

For other physically interesting choices of parameters the aim is to obtain information 
about the size of the fluctuations of the phase angle a deflned by 



We have investigated the fermion determinant in four series of points in the parameter 
space. Each series is specified by the values of G^, G^, K and the quartic scalar self- 
coupling A. The parameters are summarized in table 1. Series D has a real determinant 
and serves as a check on the program. 

Within each series the scalar hopping parameter k was varied in such a way that five 
equidistant points in parameter space are obtained, which start in the symmetric phase 
near a scalar mass of mn ~ 1 and end in the phase with broken symmetry at a scalar 
mass of the same magnitude. This covers the physically interesting region near the second 
order phase transition line. At each of the 20 points in parameter space after equilibration 
the determinant was calculated 40 times, always separated by hundred trajectories from 
each other. The results of the calculation are displayed in table 2. It displays the average 



detQ = I detQl e 



(9) 
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values of the real and the imaginary parts of det together with the average value a of 
the phase and its standard deviation 

^ 40 

^«=4oE(«i-«)'- (10) 

In scries D, where the phase is known to be zero, the standard deviation of a yields an 
estimate of the precision of the numerical results. 

From the data it appears that nothing particular happens to the determinant when 
the vicinity of the phase transition is passed. The absolute values of both the real and 
the imaginary parts increase continuously from the symmetric to the broken phase. The 
average values of a are always statistically consistent with zero, as they should. 

Most important is the observation that in all points the standard deviation of a is 
very small. Its value in series A, where it is largest, is near 2 ■ 10^^. In this series the 
Yukawa couplings differ most in their absolute values. This represents the case that is 
farthest from = iG^' where a vanishes. 

4 Conclusion 

The investigation of the complex fermion determinant in the U(1)l ® U(1)r symmetric 
Yukawa model with mirror fermions has shown that the phase of the determinant is 
fluctuating only very little, while the modulus of the determinant varies over many orders 
of magnitude in the physically interesting region. Thus a simulation of the model without 
further technical doubling of the number of fermions appears feasible. This could be 
achieved by a Hybrid Classical-Langevin algorithm based on the effective action S{(f)) — 

1/2 In det gg+. 
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Table 1 

Parameters of the four series of points at which the fermion determinant has been mea- 
sured. 



label 


A 




Gx 


K 


A 


1.0 


1.0 


0.0 


0.10 


B 


1.0 


0.3 


0.0 


0.10 


C 


10-4 


0.1 


- 0.2 


0.13 


D 


1.0 


0.3 


- 0.3 


0.10 



Table 2 

Results of the numerical calculation of the complex fermion determinant. 























label 


/V 


Re (let Q 


111! (let Q 


(\ 




A 


(\ 






0.070 


6.90 


•10^ 


2.68 


10^ 


3.89- 


10- 


-4 


2.04- 


10- 


-3 




0.079 


2.20 


•108 


-2.48 


10^ 


2.54- 


10- 


-5 


2.67- 


10- 


-3 


A 


0.088 


3.34 


•lO^i 


7.07 


10^ 


-6.17- 


10- 


-4 


2.23- 


10- 


-3 




0.097 


9.65 


•10^=^ 


-6.46 


10^° 


2.95- 


10- 


-4 


2.16- 


10- 


-3 




0.105 


7.70 


•10i« 


-2.14 


10^5 


-4.11- 


10- 


-4 


2.16- 


10- 


-3 




0.137 


1.67 


• 10^ 


4.51 


10-3 


1.15- 


10- 


-6 


1.90- 


10- 


-5 




0.145 


2.03 


• 10=^ 


5.28 


10-4 


5.55- 


10- 


-7 


2.24- 


10- 


-5 


B 


0.153 


2.64 


•10^ 


-1.35 


10-2 


-3.55- 


10- 


-6 


2.18- 


10- 


-5 




0.160 


3.09 


•10^ 


-2.10 


10-2 


-1.61- 


10- 


-6 


2.52 • 


10- 


-5 




0.168 


5.46 


•10^ 


-1.98 


10-2 


-2.34- 


10- 


-6 


2.14- 


10- 


-5 




0.077 


5.72 


.1035 


-5.38 


10=^1 


1.54- 


10- 


-5 


1.70- 


10- 


-4 




0.081 


5.77 


• 10=^^ 


-1.81 


1033 


-5.58- 


10- 


-5 


2.52- 


10- 


-4 


C 


0.086 


9.24 


.1038 


4.56 


1035 


-2.20- 


10- 


-5 


2.45 • 


10- 


-4 




0.091 


3.95 


•10^2 


-2.30 


1038 


-1.16- 


10- 


-5 


2.57- 


10- 


-4 




0.096 


2.20 


• 10^^ 


-8.98 


1060 


3.42- 


10- 


-5 


1.53- 


10- 


-4 




0.104 


3.90 


.1084 


-1.39 


10^° 


1.81- 


10- 


-16 


5.50- 


10- 


-15 




0.109 


1.10 


.1085 


-2.73 


10^0 


-7.16- 


10- 


-16 


6.53- 


10- 


-15 


D 


0.114 


3.38 


.1086 


2.52 


10^2 


9.54- 


10- 


-16 


8.06- 


10- 


-15 




0.119 


2.55 


.1088 


-1.02 


10^4 


-1.68- 


10- 


-16 


7.65- 


10- 


-15 




0.124 


1.27 


.1089 


2.44 


10^5 


2.57- 


10- 


-15 


7.42- 


10- 


-15 
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